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SUM'IARY 


A new technique was developed to systematically expand the aero the rrmodynamic 
envelope of the Space Shuttle Thermal Protection System (TPS). The tecarnique re- 
quired transient flight test maneuvers which were performed on the seccmc, fourth, 
and fifth Shuttle reentries- .Kalman filtering and parameter estimation: ^vere used 
for the reduction of embedded thermocouple data to obtain best estimates of aero- 
thermal parameters. Difficulties in reducing the data were overcome or irinimized. 
Thermal parameters were estimated to minimize uncertainties, and heating rate para- 
meters were estimated to corielate with angle of attack, sideslip, defle^rriori angle, 
and Reynolds number changes- Heating trends from the maneuvers allow fcrr rapid and 
safe envelope expansion needed for future missions, except for some loc^H areas. 


INTRODUCTION 


Because of the lifting capability of the Space Shuttle orbiter, its ranging 
capability and the aerodynamic heating to its TPS can vary significantly ~with atti- 
tude and in turn with the reentry trajectory. Flight simulators for mcsr airplanes 
today have fairly standardized equations 'of motion in terms of linearizes: stability 
and control derivatives for example. No such capability existed for aerrjcynamic 
heating- Development of a standardized procedure on flight simulators wss and is 
needed for manned reentry vehicles. 

Most heat transfer data from wind tunnel tests for the orbiter were fairly 
standardized- The ratio of film transfer coefficient to a reference sts^mation 
coefficient was tabulated as a function of angle of attack, sideslip, deflection 
angle, and Reynolds number. The wind tunnel data must be scaled to flight con- 
ditions, however, especially when the flow in the wind tunnel was transf rional . Var- 
ious theories to accomplish this were often buried in large programs whi^zni primar- 
ily output temperature time histories and were not appropriate for flimct sim-^ 
ulators or for mission planning. A simplified method (Ref. 1—2) was used for mis- 
sion planning and was adapted and modified for flight simulators (Ref. 3—5). A 
one-dimensional thermal model was used to improve accuracy for bondline nemveratures , 
and simplified heating ratios were moc ified to a tabulated form similar r-o “he wind 
tunnel ratios, or scaled ratios could be used. An added advantage of th:t_5 - jjproach 
was flexibility for updating from flight test data. 

A systems approach was used to develop the new technique for aerotz.termcdynamic 
envelope expansion of the orbiter for operational missions at Vandenbur-r Air Force 
Base. A diagram of the approach is summarized by Figure 1. The systems approach 
essentially addresses operational needs for the life of the vehicle an<f Tiot just 
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needs for the next flight. Once the vehicle design was frozen, wind tunnel data 
and thecTy established a simulator data base for heating ratios. Tabulated heatirg 
ratios and a one-dimensional thermal model for several control points (Ref- 6) 
were programmed on the simulator at the Air Force Flight Test Center (AFFTCl to 
obtain both surface and bondline temperatures. Using the thermal model also alioied 
calculation of temperatures which would be measured by existing thermocouples in the 
TPS. Thus, not only were temperatures on future operational flights predicted, bit 
thermocouple response during transient flight test maneuveijs was also simulated - 
The data reduction technique for heating estimation, referred to as HEATEST, was 
developed to estimate heating parameters in a manner similar to estimation of stai- 
ility and control derivatives during aerodynamic flight test maneuvers. Pu.<nover- 
pullup (POPU) maneuvers and flap maneuvers caused sufficient thermocouple response 
•to allow aero thermodynamic envelope expansion. The POPU allowed angle of arrack 
envelope expansion (heating as a function of angle of attack) needed for mor*e crota 
range. Flap maneuvers allowed axial center of gravity expansion (heating as a 
furction of elevon and flap deflection.). Lateral center* of gravity expansion (orat- 
ing as a function of sideslip) could not be accomplished by maneuvers. By sddresj— 
ing the overall need for flexible and quick updating of the simulator data base, 
almost identical equations for the thermal model and heating rates were used in Lie 
data reduction program. 

The data reduction program actually became a data correlation program wdien 
thermocouple data from flight test maneuvers were input. Maneuvers were designed 
to vary angle of attack, for example, while other variables were nearly consxant. 

The heating rate at a given angle of attack was assumed to be the same, and thus zj: 
correlate. Otherwise, any hysteresis during a manuever would be due to error in ts& 
thermal model. Uncertainty was thus decreased by identifying this error by estine- 
ting thermal model parameters. Use of a Kalman filter further minimized otber un- 
certainties due to modeling. Details of the technique are given in References 
and 7 - 

Lessons have been learned by applying this new technique to thermocouple dati 
from flight test maneuvers on the second and fourth Space Transportation Systems 
flights (STS-2 and STS-A). Another maneuver is also available from STS-5. Heal- 
gas effects, internal radiative and convective cooling, and late transition, are 
among heating trends which have been identified from thermocouple data, but cnly 
trends from transient flight test maneuvers will be emphasized in this paper. 


THERMOCOUPLE INSTRUMENTATION 


The excellent thermocouple instrumentation embedded in the orbiter TPS vas 
designed for typical reentry profiles, and not for transient flight test maraeuvers- 
For future vehicles, considerations for flight test maneuvers should be emp h.as izeii. 

Surface thermocouples were installed in the TPS and covered with a thin: coat- 
ing of thickness AX^. This coating was applied by hand and according to wei^t, nir 
thickness. For a transient maneuver, error in the coating thickness can cause larje 
uncertainty in the heating rate. Error in the specific heat and conductivirj of Lie 
coating could also cause more uncertainty. Therefore, an effective coating Lhickress 
(equivalent to the surface thermocouple depth and including any significant isire leac 
capacity) was estimated during flight test maneuvers by the REATEST program. 

The coating thickness of High-temperature Reusable Surface Insulation (3RSI) 
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was expected to be between 10 to 15 mils (0.025 to O.OA cm). Estimates for coating 
thickness on HRSI have varied between a maximum of 21 mils (O.Oo cm) at an outboar_ 
eleven location (V07T9730) to 15 mils at a lower surface location (V09T9527). 

Nocex felt Flexible Resuable Surface Insulation (FRSI ) coatin’ was expected tibs 
around 7 mils (0.02 cm). A location on the side of the Orbital Maneuvering System 
(OMS) pods were of primary concern because of discoloration of the FRSl during Che 
first flight CSTS-1). Because of the discoloration, the FRSI was apparently coatee 
again, increasing the coating thickness. An estimate j-or the FR_^ coa^-ng during, 
the Mach 20 POPU on 5TS-2 was 20 mils (0.05 cm) at one location (/07T99i6). However, 
a time skew of 3 seconds was also necessary for data correlation. 


Time skews were identified as a serious problem based on estimates crom sim 
ulated thermocouple data in Ref. 3. If thermocouple samples incorrectly led samples 
of the angle of attack by only one second, for example, estimates for tr.e HRSI coat- 
ing thickness became negative and physically unrealistic. If the thermocouple lagged 
by a second, then the estimates for coating thickness increaseo. Since the thermo- 
couple and angle of attack were recorded on different recorders vith no common cloc.c, 
there was and ’s a concern over time skews. The actual sample tine was also unkn^n 
within the sampling rate, which was once per second for thermocouples and angle of 
attack . 

Thermocouples were not calibrated before the firsts flights . Most of tne error 
associated with a calibration would probably only be in the form of a bias since 
calibration curves for thermocouples are well known. Most bias errors could be 
checked at ambient conditions inside the hangar within the data recorder resoluticn, 
which for the orbiter was an eight-bit word. This resolution was the primary noise 
source for the reduction of orbiter thermocouple data during transient maneuvers. 
Each thermocouple was scaled according to the anticipated maximum temperature at its 
location to minimize the resolution error. Calibration curves were then approxima..ed 
*by polynomials and an additional small error was Intioduced. 

It is suggested that thermocouple installation on future vehicles with a low- 
conductivity thermal protection system be similar to the orbiter with the following 
improvements. A pressure transducer, surface thermocouple, and condline thermocouple 
should be at the same location to enhance utilization of all measurements. A step 
input to the installed thermocouples with a known heat source should be used to 
verify the thermal model (at least at ambient conditions). The timing of the step 
input relative to the thermocouple samples is crucial to accurate estimation of 
the effective coating thickness. An accurate calibration curve could be practically 
used if all thermocouples of the same type are scaled identically and higher data 
resolution is used. Raw-data reduction would be simplified at the expense of more 
data storage capability. Real-time data links could possibly ofiset the additional 
data storage. In addition, fewer thermocouples may be necessary aecause of a better 
understanding of reentry heating gained from the thermocouple measurements on the 
orblt'-'t. More flexibility in installing thermocouples at critical locations, which 
may not be identified until after a first flight, could also reduce the number of 
measurements. Although a higher sample rate is needed during flignt test maneuver,, 
a lower sample rate may be sufficient for most of the reentry. 


THERMOCOUPLE DATA PvEDUCTION 


The reduction of thermocouple data using 


the HEATEST program was first success- 


ful7.y demonstrated by simulating thermocouple data on a typical lower surface 
location '.Ref. 3). Thermocouple data on HRSI during a transient manetzver in a wine 
tunnel test was also reduced (Ref. 3, 4, and 5), but difficulty in estimating ccat- 
ing thickness was encountered. Thermocouple data from STS-1 were reduced in Ref- 
erences 3 and 7, although data were lost above Mach 14 and there were no canei^vers. 
Thermocouple data from STS-2 and STS-4 have also been reduced (Ref. 4-6). The dat£ 
reduction technique and the lessons learned concerning the technique are discussed- 

« 

The heating equa^ Lons and one-dimensional model are referred to as the neatinr 
model and thermal model. These models were chosen to be nearly identical to the 
simulator equations. The ratio (denoted by a bar) of the heating rate (q) or film 
transfer coefficient (h) to an appropriate reference condition was assumed to be a 
linear function of the form 


^ ^log RE RE -_iog RE^) 

^ q^e ^ 

where q is the magnitude or intercept at the reference conditions specified by the 
zero suSscript on each variable. The subscripts on the heating ratio (q) represent 
partial derivatives with respect to each variable. The variables are angle of attaik 
(a), sideslip (S), logarithm of the freestream Reynolds number (log RE), eleven de- 
flection angle (6e) , and flap deflections • Second-order terms wi.th second-deriv 

ative parameters vere also added somtimes to account for nonlinearity. The heating 
model in the simulator was analogous to the linear equation when the heating param- 
eters q 9 q f qns q.j^Q j^9 ^ibf functions of the appropriate variables- 

For data^reduct ion dur?ng maneuvers, these parameters were assumed to be constant ter 
short time durations, but to vary during the reentry. The heating rate was obtained 
by multiplying the ratio by the reference heating (q ). The heating rate was input 
to the one-dimensional thermal model. A typical TPS cross section for Reusable 
Surface Insulation (RSI) is shown in Figure 2. The thermal -mode I equation.^ were 
solved by the same finite-element (or finite -diffference) method as on the simulator 
and require initial conditions. 

A simplified diagram of the HEATEST algorithm is shown in Figure 3- The 150DEL 
block was identical to the simulator model except the sensitivity and covariance of 
the temperature at each node were also propagated in time. Initial conditions in 
the IC block were required for temperature (U), the sensitivity temp- 

erature to each parameter (6 ), and the covariance (P) of the temperature. The 
parameters (9 ) include bor.h^thermal and heating model parameters. The propagated 
or predicted lEemperature , sensitivity, and covariance are referred to as a priori 
values and are denoted by a minus superscript. 

If a thermocouple measurement was available, then the a priori temperature was 
compared with the measured temperature in a Kalman filter algorithm (Ref. 3, 4, 5, 
and 7) which is referred to as a KALMAN UPDATE. Depending on the Kalman filter 
tuning (i.e., the uncertainty in the measurement and models), the a priori values 
were updated to a posteriori values which were denoted by a plus superscript. T’his 
process was repeated in a TIME LOOP until the end of a given time segment.- 

A Newton-Ralphson algorithm was then used to update heating and thermal para- 
meters in the PAR/J-IETER UPDATE block to satisfy a maximum likelihood function for 
each parameter (Kef. 3, 4, 5, and 7). With the updated parameter, the TIME LOOP 


was rtpeated. This process was repeated for a fixed number o: iteraricns, 


FLIGHT DATA REDUCTION 


: r 1 c u 
r co:v 


rry 
rmoc:'. 
I- s 

lost : 
Ivity. 
ronct*: 


The HEATEST algorithm was applied to flight thermocouple lata. :>ome 
was encountered due to data loss, initial condition genera-^ lo. , and .xce^si 
puter time. 

Onboard thermocouple data were lost on STS-1 and STS--, :^.*d te~ex^ 

data were obtained. Therefore, the only flights with maneuvers ana vita th^ 
data are STS-2 and STS-5. Data from STS-5 were not availaole in time for t.. 
paper, but do not appear to change any results. Pressure meaiureme,_^s ^ 

STS-2. Therefore, predicted pressures were used for determining RS - conduce 
Telemetry data from STS-4 partially covered both FOPU maneuvers, but cause 
for initial conditions. 

Errors in initial conditions were determined to be a source of error in nea.-.g 
parameter estimates based on studies with simulated thermocouple data y-ef. ^ 

STS-2, initial conditions were generated by reducing all tr.errocoup.e data 
interface' to just prior to the maneuvers at approximately !lacn _ . -ni.-^ 

on-board dynamic pressure and velocity were used to calculate - 

pitch jets fired at low dynamic pressure, a large spike was c..used ‘ - 

pressure and thus in density. HEATEST estimated large angle or attar^ 
since no input variables except angle of attack were changing. If ^ro'-a- 

were input and the heating allowed to correlate with this var...ble, ...^n “ 

should identify the error. Time segment size was increaseo insteaa.^ I. " 

conditions for STS-2 maneuvers should therefore be accurate?, out excess* 
time was required. . * 

Several methods to decrease the computer tiiie \ J^°"Ul-AdarT 

time was used to propagate the covariance. First : 

condition was assumed only for the covariance at a few noces 

Another approximate initial condition generation procedure ^ i;. 

analogy with an ^erically determined time constant (Ref. 4 o). Y. 

efficient libr routines were replaced by new f’’^%=f‘"^;®“L^ 57 -VD 7 T 953 irwitv 
Analysis of th ach-20 POPU at a lower surface location (709-^.27 4 D7T9^3 ; 

14 requited 440 seconds on a Cyber 74. The nen "“j"" 

to 44 seconds. Using the adiabatic «all assumption would 

would not alio. th. Ss. of .11 lo-depth thetmocoup l.s . f f:r“;,tla 

tlons w.te generated 50 seconds prior to the =d. 

conditions for the POPU. The behavior of the Kalman filter a. ea 
was demonstrated as shown in Figure 4. 


TRENDS FROM FLIGHT TEST MANEUVERS 


Heating and thermal parameters were estimatea during ^r.uve . s 
effective coating thickness have already been discusse . „ea -,.g 
maneuvers were obtained at a few locations and 

ative of the lower surface, the eleven control surface and ).u. pods .... .? 

surf ace . 



Lower Surface 


because of the data loss on STS-4, a Icwer-surface location (V09T9527-V07t9531 - 
with, five embedded thermocouples was investigated. Tine histories of the surface 
thermocouple and angle of attack of the Mach 20 POPU during STS-2 are shewn in 
Firrre 5. Time histories of the Mach 12 and 8 POPU during STS-4 are shown in Fig- 
ures 6 and 7. Note that the pullup portion of the Mach 12 maneuver was unusually Icm; 
and chat telemetry data started during the pushover. Note theC the Mac]i & m.aneuver 
was short duration. 

?.£sults for the Mach 2C maneuver agreed well with wind tunnel data as shevn in 
Figure 8. The trend in angle of attack agrees well and demonstrates that the 
linear assumption in angle of attack was valid. 

The angle of attack heating derivative was lower at Mach 12. This coula he a 
Mach effect or initial condition error caused by the data loss and transition to 
fully turbulent flow above 44 degrees angle of attack. As seen in Figure 6, the 
lennerature was higher than the simulated temperature based on laminar heating rale^ 
In addition, as suggested by Hertzler (Ref. 6), the axial accelerometer nicasure- 
meni (axial drag) increased discretely for angles of attack above 44 degrees, typi- 
cal of flow transition, and decreased back to laminar flow below 44 degrees. Thus, 
sensitivity of transition to angle of attack was demonstrated. 

The flow was fully turbulent during the Mach 8 POPU. The higher heating nag- 
niende and larger uncertainty bound are shown in Figure 7. The higher uncertainty 
was due to the short duration maneuver and perhaps to heating changes with Reynolds 
nun±er. A Reynolds number derivative could not be estimated even when the angle of 
attack derivative was fixed. An alternate, procedure used sequential five second 
time segments- The trend in Reynolds number, assuming the angle of attack trend wa^ 
correct, is shown in figu.e 9. 


Zlevon Control Surface 


Beating trends were sucessfully eclimated at a location near the tip of tb* 
outbeard elevon (V07T9730) for numerous maneuvers during STS-2. Time histories are 
shown for the Mach 21 flap maneuver and Mach 20 POPU in Figure 10. Because there 
was no thermocouple response, the flap maneuver at Mach 16 is not shown. The >:ach 
12 flap maneuver is shovm in Figure 11. 

Nonlinear beating trends in elevon deflection angle were evident as shown in 
Figure 12 for rhe Mach 21 maneuver. A second-order polynomial in elevon deflection 
was ’^serd. The cause of th€r increased heating above five degrees elevon deflection 
is unknown but could be due to a localized flow phenomenon or to transition to 
turbulent flow. There is a possibility that a local separation bubble would raise 
the hypersonic buffet reported on STS— 4 when the elevon schedule was five degrees. 
Beca-cse there ves no thermocouple response during the Mach 16 maneuver, the eleven 
heating derivative was zero. This result implies that the heating estiiiates at 
Mach 21 are more nonlinear. 

P^asults fron the Mach 12 flap maneuver for fully turbulent flow on the elevon 
are also shown in Figure 12- The data correlated poorly for negative deflection 
angles. A localized phenomenon or transitional flow at small deflections could 
cause this behavior. Fcr envelope expansion, the tre:.ds at the larger deflection 


angles were of primary concern- The lower e^'cvon heating allowed for a traLic— jlc 
with the flap control surface. 

Hearing on the llap was lower than expected for a gi'ren de fleet icr* angle, bu- 
because of a shift in the basic pitching nsoment, larger ilap or elevc:c deflection. 5 
were required. ihe heating on the flap changed dramatic.= lly with fla^ defle^cic., 
angle of attack, end elevon deflection during STS-2 nanervers. There 'was al^o ci 
Reynolds number trend which made data correlation very c^fficult (Rc-» V). 


Orbital Maneuvering System ?cc.s 


The primary concern for angle of attack envelope expansion was the upper sur_ecc, 
especially the OMS pods. Time histories for the Mach 20 ?OPU during STS-2 are sr.-n-wn 
in Figure 13. A large unexpected response occurred on t r e side of the OMS pod 
(V07T9976). Although a three-second time s'<ew was necessary, the results from param- 
eter estimates correlate well as shown in Figure 14. The flow impingement on the 
pod started at 37 degrees angle of attack instead of at 30 degrees as in the wine: 
tunnel. In addition, estimates in Ref. 6 from STS-4 indicate less decadence on 
Reynolds number than in the wind tunnel. Visual inspection of FRSI discoloration 
as shown in Fig. 15 indicated a different pattern than expected based on predict iims 
from wind tunnel data. 

Investigation of transient maneuvers in the wind tunnel in Ref. 1C to estimate 
coating thickness demonstrated a similarity with the OMS pod heating.^ A three-second 
lag in the thermocouple response was found. A fheoretic^il investigation in Ref, i-l 
of the lower heating magnitude reported in Ref. 3-5 confirmed the sensitivity of 
the heating to a discontinuity in wall temperature (cr ncmisothermal wall caused p' 
an interface between different materials). The stainless-steel leading edge and .-RST 
test article on a flat plate in the wind tunnel had a step increase in wall temp- ^ 
erature. On the CMS pod, there was a step decrease at low- temperature rsi ancf 

FRSI interface as showTi in Fig. 15. Therefore, the taree— second tine skew and semte 
of the increased heating on the OMS pods mar be attributed to the nonlsothemal wa.il. 

Numerous locations on the Orbiter have an interface between different matoria-s 
and probably have the same nonisothermal-wall problem. Coe such local area would oe 
the nose cap and HRSI interface where tiles have sluiced due to increased heating 
as shown in Fig. 16. According to Ref. 11, such an irterxaca near the leading edge 
would have a rapid recovery. The nose cap surface temperature peaks around Mach -3 
corresponding to a POPU on STS-5. The nonisothermal wall effect woulc be largest at 
Mach 18. This effect should be considered when heating c-ata is correlated and aL*.o 
for future designs. 
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Figure 1.- Systems approach to envelope expansion. 
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Figure 3.- Simplified HEAIEST flow dlagrsm. 
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Figure 7.- Lower-surface-plug STS-4 flight thermocouple 
data for turbulent flow (Mach 8 pullup-pushovcr 
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Figure 8.- Heating estimates for lower surface plug 
from STS-2/STS-4 flight thermocouple data. 
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Figure 9.- Heating comparison for STS-A thermocouple data 
using both stagnation and Eckert turbulent reference 
heating to show Reynolds number trends. 



Figure 10.- Outboard eleven STS-2 flight thermocouple 
data (Mach 21 flap maneuver and Mach 20 pushover- 
pullup maneuver) . 
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Figijre 11.- Outboard eleven STS-2 flight thermocouple 
data (Mach 12 flap maneuver). 
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Figure 12.- Heating estimates for outboard eleven 
frc« STS-2 flight thermocouple data. 
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Figure 13.- OMS pod STS-2 flight thermocouple data 
(Mach 20 pushover-pullup naneuver) . 
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Figure 14.- Heating estimates for OMS pod 
from STS- 2 flight thermocouple data 
(Mach 20 pushover-pullup maneuver) . 















